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Abstract 

We consider a physical system with a couphng to bosonic reservoirs via a quan- 
tum stochastic differential equation. We study the limit of this model as the cou- 
pling strength tends to infinity. We show that in this limit the solution to the 
quantum stochastic differential equation converges strongly to the solution of a 
limit quantum stochastic differential equation. In the limiting dynamics the ex- 
cited states are removed and the ground states couple directly to the reservoirs. 



1 Introduction 



It is a frequent occurence in physics to have a system that spends a very limited amount 
of time in its excited states. This is, for instance, the case if the system is strongly 
coupled to a low temperature environment (e.g. the optical field). The strong coupling 
ensures that excitations above the ground levels of the system quickly dissipate into its 
environment. It is therefore reasonable to ask for a model in which the excited states are 
eliminated from the description. That is, we would like to have a description that only 
involves the ground states of a system and its environment. The procedure for going 
from the full model to the reduced model is called adiabatic elimination. 

We study adiabatic elimination in the context of quantum stochastic models |15| which 
arise by taking a weak coupling limit of QED (quantum electrodynamics) models pQ 
[T3l [5], and are widely applicable to systems studied in quantum optics. Specifically, 
quantum stochastic models are the starting point for deriving master equations, filter- 
ing equations, and input-output relations. In the quantum optics community adiabatic 
elimination is a common technique, used, for instance, in atomic systems [22l O |6l [11] 
and in cavity QED models [HI [1^ as well as in more recent work on quantum feedback 
[71 [9l [24] . Rigorous results have been demonstrated for adiabatic elimination outside of 
the quantum stochastic models we consider [Tni[ll[II]. At present, however, apart from 
the work [14] on the elimination of a leaky cavity (using a Dyson series expansion to prove 
weak convergence), no rigorous results have been obtained on adiabatic elimination in 
the context of the quantum stochastic models introduced by Hudson and Parthasarathy 

m. 
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We start by considering a family, indexed by a parameter k, of quantum stochastic 
differential equations (QSDE's). The parameter k can be interpreted as the coupling 
strength between the system and its environment. The environment is modelled by a 
collection of bosonic heat baths in the vacuum representation. We assume that the 
coefficients of the QSDE are all bounded and satisfy the usual conditions guaranteeing 
a unique unitary solution |15j . We state further assumptions on the coefficients and 
show that under these assumptions the solution of the QSDE converges strongly to the 
solution of a limiting QSDE as k tends to infinity (Theorem 12. ip . The limiting QSDE 
represents the adiabatically eliminated time evolution of the system. 

The heart of the proof is a technique introduced by T.G. Kurtz [T71 that enables the 
application of the Trotter-Kato Theorem [5T] . This allows us to prove strong convergence 
of the unitaries using convergence of generators of semigroups rather than convergence 
of a Dyson series expansion. Convergence is first shown on the vacuum vector of the 
bosonic reservoirs. We then extend this result to any possible vector in the Hilbert space 
of the reservoirs by sandwiching the unitaries with Weyl operators and using a density 
argument. 

The remainder of this article is organized as follows. In Section[2]we introduce the system 
coupled to n bosonic reservoirs in the vacuum representation. We state assumptions on 
the coefficients of the QSDE and present the main convergence theorem. In Section 
[3] we discuss four applications of the theorem in the context of examples from atomic 
physics and cavity QED. Section 3] presents the proof of the main convergence theorem. 
In Section [5] we discuss our results. 



2 The main result 



Let 7i be a Hilbert space and let n be an element of N. Let T be the symmetric Fock 
space over C" (g) L^{R+) = L^{R+; C"), i.e. 

OO 

J- = c e l2(r+; c")®-'". 

m— 1 

Physically, the Hilbert space TL®!F describes a system Ti. coupled to n bosonic reservoirs 
(e.g. n decay channels in the quantized electromagnetic field). For / g L^(M+;C"), we 
define the exponential vector e{f) in J-' by 

e(/) = l®0:^. 

m— 1 ^ 

Moreover, we define the coherent vector ■n{f) to be the exponential vector e(/) normal- 
ized to unity, i.e. 7r(/) — exp(— i||/|p)e(/). The vacuum vector is defined to be the 
exponential vector <& ~ e(0) = 1 © © . . .. The expectation with respect to the vacuum 
vector is denoted by 0, i.e. </> is a map from B{!F) (the bounded operators on J^) to C, 
given by (l){W) = ($, W^) for aU W G B{T). 
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The interaction between the system and the bosonic reservoirs is modelled by a quantum 
stochastic differential equation (QSDE) in the sense of Hudson and Parthasarathy [TS] 
of the form 

dC/f-) = I - dAj^- + if - Lf^^sg^dAi + K^'^^dtY^'^ , (1) 

(k) 

where Uq = I. We consistently use the convention that repeated indices that are not 
within parentheses are being summed (i and j run through {l,...,n}). The Hilbert 
space adjoint is denoted by a dagger We have indexed the equation with a positive 

(k) 

number fc, and in the following we will be interested in the behaviour of as k tends 
to infinity. We assume that the following conditions on the coefficients of the QSDE are 
satisfied. 

Assumption 1: For each A: > 0, the coefficients K''-'^\ s['^^ and hf^ e {1, . . . , n}) 
of the quantum stochastic differential equation ^ are bounded operators on Ti.. Fur- 
thermore, for each fc > 0, the coefficients satisfy the following relations 

Hudson and Parthasarathy [15j show that under Assumption [H the quantum stochastic 
differential equation ^ has a unique unitary solution , and, the adjoint satisfies 
the adjoint of Eq. 

Assumption 2: There exist bounded operators Y, A, B, Fi, Gi and Wij (independent of 
fc) on H such that 

K(k) ^pY + kA + B, Lf^ = kF, + , 4^"^ = ' 
for all i, j G {1, . . . , n}. 

We define Pq as the orthogonal projection onto Ker(y). Let Pi = / — Po be its comple- 
ment in Ti. We use the following notation Tio = PqH and Tii — PiH. Physically, one 
should think of Tio as the ground states and of Tii as the excited states of the system. 

Assumption 3: There exists a bounded operator Y^"^ on Ti. such that PiY^^ — Y^^Pi 
and 

yyf iPiZPo = PiZPo, PoXPiY^^Y = PoXPi, (2) 

where Z = A,F^W^,, (j e {l,...,n}) and X - A, P, P„ G„ W,, , ClWij , P.Ff ^P, , 
F,Y,-^A,F,Y,-^F^Wij,AYf^A,AYf^F„AYf^F^Wi,, e {l,...,n}). Moreover, for 
alH, j e {1, . . . , n} the following products are zero 

PoyPi = PoAPo = P,Po - PoiSa + F,Y^-^F^)Wi,Pi - 0. 
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Note that the existence of ^ satisfying the assumptions in Eq. ^ is immediate if Y 
maps Til surjectively onto Tii and is therefore invertible on Tii. 

Definition 1: Suppose Assumption [2] and [3] hold. We define for all i,j £ {1, . . . ,n} the 
following bounded operators on Ti 



K^Po{B~ AY^^A) Po, U - (G. - F^Y^^A) Pq, 



Assumption 4: For all i, j e {1, . . . , n} the following products are zero 

PiL, = PiS,j = 0. 

Lemma 1: Suppose that Assuinption\^\^\^and\^hold. The operators in DeHnition\^ 
satisfy 



Proof. By Assumptions[I]and[3we have ifC^) = -hf^^hf', K^^'^ = k^Y+kA+B 
and i^'^'' = kFi + Gi for all k >Q. Moreover, P^Pq = 0, by Assumption [3] Combining 
these results leads to 



We then use YY^ ^^Po ~ AP^ from Assumption [3] and Li from Definition [T] to derive 



= -Po{g\ - a^y-'^f}){g, - p.yf M)Po 

= Po(P + Pt)Po _ PoAty-it(^ + ^t)p^ 

- Po{A + At)y-iAPo + PoAt(y-it + Y^^)APo 
= Po(P + pt)Po _ PoAFf MPo - Po^lty-it^tp^ 
= Po(i^ + /v^)Po. 



Combining this with —F^Fi = Y + Y^ from above, 




p/Pi = y + 

PiP;G,Po-Pi(A + At)Po 
PoG^G,Po = Po(P + pt)Po. 



(3) 



By Definition [U 





PoWlWijPo = Po6, 
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Then use Pq [Sa + F^Ff ^f/J Wi-jPi = from Assumption [3] and PiSi-jPo = from 
Asumption [4] to derive 

SaSl = Po (<5„ + F,Y^-'Fl) WniW^i (Smj + F^Y-^^ f]^ Pq 
= Po (<5.„ + F,Y^-'Fl) (5„, + F„r-i^F/) Po = %-Po. 

□ 

The operators given by Definition [T] are the coefficients of a QSDE on the Hilbert space 

dUt = i^[s,j-S^jPo)dA't^ +L,dAl^ -Lls^jdAl +Kdt^Ut, Uo = L (4) 

Lemma [T] imphes that under Assumptions [1] [51 [3] and 31 Eq. ^ has a unique unitary 
solution on H [T5], and, the adjoint u} satisfies the adjoint of Eq. Moreover, Ut 
maps Hq to Ho- Note that UtPi — Pi. 

Theorem 2.1: Suppose Assumption [210 El and [3 hoW. Let u[^^ he the unique unitary- 
solution to Eq. Ilj. Let Ut be the unique unitary solution to Eq. ([4|) where the coefficients 
are given by Deenition\^ ThenU^'^^Po converges strongly to UtPo, i.e. 

hm [//''V = Ut^, VV' e Ho ® J^. 

k — ^oo 

We prove Theorem 12. fl in Section [H 



3 Examples 



We use the fohowing definitions in the first two examples below. Let (|e), \g)) be an 
orthogonal basis of C^. Define the raising and lowering operators in this basis as 

fO l\ fO 

"+ = [0 oj' "-^[1 

Define the Pauli operators 

ax = 0"+ + o"— , cTy = — 'i(T-|_ + icr_, (T^ = (T+(t_ — (7_(t_)_, 
and define the projectors 

Pe=(T+(T_, Pg^a-a+. 
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Example 1: (A two- level atom driven by a laser) The Hilbert space for a two-level 
atom is 7i = C^, with |e) the excited state, and \g) the ground state. Define the detuning 
A € K, the decay rate 7 > and the complex amplitude a G C. The QSDE for this 
system in the electric dipole and rotating wave approximations is [5] 

dUi'^^ — fc^/7cr_ dAj — ky/^a+dAt — ikacr+dt — ikda-dt 

- ^a+a_dt - ik^Aa+a_dt^ui''\ U^''^ = I. 

Define the operators Y, A, B. F, G, W as 

Y = (— iA — 7/2)(T-|_cr_, A = —iaa^ — iaa-, B = 0, 
F = y/ja-, G = 0, W^I. 

This satisfies Assumptions [T] and [U and Pq = Pg. We take Y^^ = —{iA + j/2)~^a^a^, 
and Assumption [3] holds by inspection. Definition [T] leads to the following coefficients 



r-i . Q^\/7 iA — 7/2 

^""zA + 7/2^^' ^ = "'zA + 7/2^^' ^ " TaT^^'' 

Note that PiL = PiS — satisfying Assumption^ Theorem 12 . II then shows that u[^'^ Pq 
converges strongly to UtPo, given by 

li. 

zA + 7/2 



dUt = . . ,„ |-7rfAt - iay/jdAj + ia^dAt - \a\^dt^ Ut, Uq = /. 



In the case that 7 = the two level atom decouples from the field. In this case we may 
explicitly calculate the ground state evolution as 

p^^~^{kaa^+k&a^+k Aa+a^)tp^ ^ _ {ft COs{knt) + ik A Sm{knt)) , 



with = ^ A^k'^ + 4|ap. For fc — > cx3 this expression limits to e*'"' which is the 

lap 

solution to our eliminated differential equation dUt = i-^Utdt^ Uq = I. 

Example 2: (Alkali atom) Now consider a system with Hilbert space H — C"^ <E) ■ 
Physically, the system represents an alkali atom with no nuclear spin coupled to a driving 
field on the S1/2 P1/2 transition. We have four orthogonal states in this system 
corresponding to the atomic excited and ground states with angular momentum ruz = i 5 
along the z-axis. We define a detuning A g R, a decay rate 7 > and a magnetic field 
-Bi G R, id x,y,z. The system may emit into n = 3 independent dipole modes, Al, 
where the modes are labelled by z G {x,y,z}. The QSDE for this system in the dipole 
and rotating wave approximations is [2] , 

dU^''^ = jfc^cr- (g) (TidAf - ky/j(7+ (g) (JidAi - ^y^-Pe <g Idt 
- i {k^APe g) / + / g) B,ai) dt \ f/f ^ , C/^''^ = /. 
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Defining the operators Y, A, B, Fi, Gi, Wij as 



Y = I -iA Pe<»I, A = 0, B ^-il® B.Ui 



2 

Fi = ^/7o•- Gi ^ 0, Wi-j = (5^ , 

satisfies Assumptions [D and [21 and = Pg ® I. We take ^ = -(iA + '^y^P^ /, 
and Assumption [3] holds by inspection. Define the eliminated coefficients as 

K = -iPg (g) BiCTj, Li = 0, Sij ^ Pg(E) / - ^ 3^ '^i'^j^ ■ 

This satisfies Assumption E) Theorem 12 . f I then shows that U^.^^ Pq converges strongly to 
UtPa, given by 

dUt = Pg®\ ~iB,a,dt - —L^aiajdAl^ | Ut, Uo = I. 



In the following two examples we make use of a truncated harmonic oscillator. We have 
truncated the oscillator to satisfy the boundedness condition of Assumption [1] in the 
following two examples. Let N be an element in N such that N > 2. The Hilbert space 
of the oscillator is C^. We choose an orthonormal basis (|0), . . . , |iV — 1)) in C^. The 
annihilation operator b : — > is given by 

6|n) = 1), n G {1, . . . , TV - 1}, 

and 5|0) = 0. The creation operator is defined to be the adjoint b^ . 

Example 3: (Gough and Van Handel [14] ) Let f) be a Hilbert space. We define 
7i = f) C^. The Hilbert space t) describes a system inside a cavity. We model the 
cavity as a truncated oscillator C^. Let Eij, i,j E {0, f} be bounded operators on f) 
such that eJj ~ Eji and ||i?ii|| < ^. Consider the following QSDE 

= l^^kbdAl - ^kb^At - ^b%dt - iH^'^^dt^ ijf \ C/($^'^ = /. (5) 

Here 7 is a real parameter and H^^'^ is given by 

^ k^Eiib^b + kEwb'^ + kEoib + £^00- 
Define operators Y, A, B, F, G, W as 

Y = [~iEii - I) 5t6, A = -i {Ewb^ + Eoib) , B = -lEoo, 
F ^ ^b, G = 0, W = L 

This satisfies Assumptions [T] and [2] and Pq = A ® Since < ^, the inverse 

(iEii + ^ exists. Let N^^ : Hi Hi be the inverse of the restriction of b^b to Hi. 
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Taking Yj^ ^ = — (iEu + -j) ^ ^Pi satisfies Assumption [3l Definition [T] leads to the 
following coefficients 



(6) 



lEii + ^ 

—iy/^l iEii — ^ 

L = — — -^EioPo, S = — — -^Po- 
lEii + \ lEii + ^ 

These coefficients satisfy Assumption 21 Theorem 12.11 then shows that U^^^ Pq converges 
strongly to C/tPoj where Ut is given by 

dUt = {(5 - Po)dKt + LdA\ - L'^SdAt + Kdt^ Ut, Uq = /. 

Remark 1: Note that we consider a truncated oscillator, where [2] treats the full 
oscillator, and that we prove our result strongly, whereas 14J proves a weak limit. The 
convergence of the Heisenberg dynamics follows immediately from our strong result. 
Apart from these points. Example [3] reproduces the result in [14j . Care must be taken 

when directly comparing the limit equations, since the results in [14j are presented in the 

(k) 

interaction picture with respect to the cavity. Under our assumptions, we define as 
the solution to 



^kbdAl - ^kbUAt - ^6tMi I f/'^) , I/q^''^ 



The unitary in the interaction picture is then given by = where t/j'^' 

is given by Eq. Note that due to Theorem 12.11 V^'^^^Pq converges strongly to VtPo, 
where Vt is given by 

dVt = ~2PodAtVt, Vo = /. 

This accounts for the sign difference between the coefficients in the equation for Ut 
presented in pj , and the coefficients in the equation for Ut given by Eq. 

Example 4: (Duan and Kimble We again consider a system inside a cavity, 
described by the Hilbert space TC = i)^ C^. The system inside the cavity is a three level 
atom, i.e. (} = C'^. Let (|e), |+), |— )) be an orthogonal basis in (). In this basis we define 

Moreover define and P± — The QSDE for a lambda system with 

one leg (+ ^ e) resonantly coupled to the cavity, under the rotating wave approximation 
in the rotating frame, is, 

dC/f ^ = J ^kbdAl - ^kbUAt ~ ^b'fbdt + 





1 


















• 









• 



















5fc2(aW6 _ a^+h^yt + k (a^T^a - a^T^a) dt C/f \ U^ 



(k) _ J 
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Here 7 is a positive real parameter and a is a complex parameter. Note that we extend 
the model from |8] to allow driving on the uncoupled leg (— <-!■ e) of the transition. 
Define operators Y,A,B, F, G, W as 

F = y^b, G = 0, W = L 

This satisfies Assumptions [T] and [2] and Pq = (|+)(+| + |->(-|) «> |0)(0|. We define 
the following subspaccs of H 

Hn = span||+) ® \n), |-) (g) |e) (g) \n - n G {1, . . . , TV - 1}, 

AT = span||e) (g) |7V- 1>|. 



7? 



Note that Hi = ^n=i and that the subspaces i?„ (n G {1, . . . , A''}) are all invariant 
under the action of Y. On the subspaces Hn, n S {1, . . . , iV — 1}, y is given by 



Y\h^ - I 






with respect to the basis (|+) (g |n), |— ) (g) |e) (g |n— 1)). Moreover, YIh^, 
The inverse is readily computed to be 

Y\-H=—A ^ 1, nG{l,...,iV-l}, 



_ 7(jV~l) 



2 




where d = ^ + ^^rt. Moreover, Y\^^^ = " ^(jv-i) • -'^'-'^ define Yj^ ^ = 

ffi^=i^l/f^ -Pi- This satisfies Assumption [3] Definition [1] leads to the following coeffi- 
cients 

^ = -^P-® |0)(0|, L=-^aL+V^^)(g|0)(0|, 5 = Po-2P_(g|0)(0|. 

These operators satisfy Assumption 01 Theorem 12 . f I then shows that U^'^^ Pq converges 
strongly to UtPo, where Ut is given by 

dUt = [{S- Po)dAt + LdA\ ~ F^SdAt + Kdt^ Ut, Uo = F 

Note that the ground state system is a two- level system on which S acts as cr^ . 



4 Proof of Theorem 12.11 



Definition 2: Suppose Assumptions [U [51 [3] and d] hold. Let B{H) and B{Hq) be the 
Banach spaces of all bounded operators on H and Hq, respectively. We define for all 
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t > and fc > 



Tt 
Tt 



(fe) 



(X) = id ® (c//xt/f ^) , X £ Bin), 

{X) = k\®(t>{ulxUtY XeB{Ho), 



where U^'^^ and Ut are given by Eqs. ^ and ([4]), respectively. 



(fe) (fe) 
Note that is intentionally skew with respect to Ut and . 



Lemma 2: For eaci fc > 0, the families of bounded linear maps Tf^'^ it > 0) and Tt [t > 
0) given by Definition El are norm continuous one-parameter contraction semigroups with 
generators 



S£'^^\X) = K^X + XK'^^^ + L\XLf\ X e B{n), 
^{X) = K'^X + XK + l\xL,, XeB{Ho), 

respectively. That is T^^'' = exp (iJf C^^) and Tt = cxp(iJf ) for aii t > 0. 



(7) 



(fe^ 

Proof. We only prove the lemma for Tj . The proof for Tt can be obtained in an 

analogous way. Since the conditional expectation id (8) (^i is norm contractive and Ut and 

(fe) 

[/( are unitary, we have 



j,(fe)(^) 



< 



u}xu, 



(fe) 



< 



u 



u 



(fe) 



\x\ 



\x\ 



for all X G B{T-C). This proves that T^^^ is a contraction for all i > 0. An application of 
the quantum Ito rule [15] , together with the fact that vacuum expectations of stochastic 
integrals vanish, shows that 



dT^^^'^ (X) = id ® ( d ( UjXU^ 



id( 



Uj (K^X + Xif + LjXLf n C/f M dt = T^ 



.(fe) 



dt. 



for all X e B{H). That is, t/ ■* = exp (t^'^'')) is a one-parameter semigroup with 
generator J^'^^^ . Furthermore, J^'^^'^ is bounded 



££^^\X) 



< 



-(fe) 



1^11 



which proves that T^'^^ is norm continuous. 



□ 



The proof of Theorem 12.11 relies heavily on the Trotter-Kato theorem [21] [16] in com- 
bination with an argument due to Kurtz [I Jj. We have taken the formulation of the 
Trotter-Kato theorem from [31 Thm 3.17, page 80], see also 10, Chapter 1, Section 6]. 
The formulation is more general than needed for the proof of Theorem 12.11 
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Theorem 4.1: Trotter-Kato Theorem Let B be a Banach space and let Bq be a 

closed subspace of B. For each k > 0, let T^^^ be a strongly continuous one-parameter 
contraction semigroup on B with generator . Moreover, let Tt be a strongly contin- 
uous one-parameter contraction semigroup on Bq with generator ^ . Let V be a core for 
I£ . The following conditions are equivalent: 



1. For aUX eV there exist X'-'''^ e Dom (^C^)) such that 



lim XC^) = X, lim ^C^') (x^'^A ^ ^{X). 

k — >oo k — *oo \ / 



2. For aUO< s < oo and all X £ Bq 



lim < sup 

k~*co I o<t<s 



T^''\x)~Tt{X) 



Proposition 1: Let t/'^' and Tt be the one-parameter semigroups on B(Ti.) and B(Ti.o) 
defined in Definition [21 respectively. We have 



lim < sup 

k^ao \ o<t<s 



Ti''\X)-Tt{X) 



0, 



for all X e B{no) and < s < oo. 



Proof. The proof follows the line of the proof of [TTl Theorem 2.2]. Lemma [5] shows 
that T^^^ = exp (t^^'"')) and Tt — exp(t^) are norm continuous, and therefore also 
strongly continuous semigroups with generators given by Eq. ([7]). This means we sat- 
isfy the assumptions of the Trotter-Kato Theorem (Thm. 14. ip with T> = B{T-Co) and 
Dom(^('=)) =B{n). 

We can write = ^oiX) + k^i{X) + k^^2{X), X e B{n), where (recall 

Assumption [5]) 

^o{X) ^ K'^X + XB + lIxG^, ^i{X) = XA + lIxF,, ^2{X) = XY. 

Let X be an element in B{'Ho) and let Xi and X2 be elements in B{'H). We define 
= X + j:Xi + ■^X2. Collecting terms with equal powers in k, we find 

J^W (xW) = (^o(^)+ =^1(^1) +=^2(^2)) + 
k{^,{X)+^2{Xi)) + 

k^^2{X)) + 



i (ifo(^i) + ^1(^2)) + ^ (=5^0(^2)) . 



Note that ^2{X) ^ as X e B{Ha) and P^Y = 0. Using the existence of \ we set 

Xi - -^i{X)Y-^Pi, 
X2 = - (ifo(^) + ^i(Xi)) Y^^Pi. 
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Using the properties of ^ in Assumption [31 we obtain 



lim (X^'^^) - _lim (^{X) + i (^o(^i) +^1(^2)) + ^^o{X2) 



fc— >oo 



A;—* 00 



The proposition then foUows from the Trotter-Kato Theorem. 



□ 



Note that for all v E Ho, we can write Utv $ = PoUtv (X) $. This leads to 



(C/f ^ - Ut)v (g) $ 



(^;,(2/-TW(Po)-T('=)(Po)^)^;) 



Here we have used that id (8)0 is a positive map, i.e. it commutes with the adjoint. Using 
Proposition [T] and noting that J^{Pq) = by Lemmas [2] and [U we see that Theorem 12. II 
holds for all vectors in Tio ® T oi the form -0 = w $. We now need to extend this to 
all il^ET-lo^T. 

Let / be an element in L^(M+;C"). Denote by ft the function / truncated at time i, 
i.e. /t(s) = /(s) if s < t and /t(s) = otherwise. Define the Weyl operator W{ft) as the 
unique solution to the following QSDE 



dWift) = { f{t),dAl^ - f{t)^dAl - lf(t)J{t),dt } W{ft), Wifo) 



(8) 



Note that W{ft) is a unitary operator from !F to !F. Moreover, it is not hard to see 
that 7r(/t) = W{ft)^, see e.g. [30]. Often we will identify a constant a e C" with the 
constant function on R"*" taking the value a (truncated at some large T > so that it is 
an element of L'^ (R+ ; C")) . 



Definition 3: Let / be an element in _L^(R+;C"). Suppose that Assumptions [1] [2l [3] 
and [4] hold and let U^''^ and Ut be given by Eqs. ([T]) and (j4|), respectively. Define 



T}''^\x) = id®^(uy^^XU, 



(kf) 



{X) = id (g) (u^^^xu^'' 



Definition 4: Let a be an element in C" and let i be an element in {1, . . . Let 

K^''\ K, L'f \ Li, S'^j ^ and Sij be the coefficients of Eqs. ([T]) and Q. Define operators 

^(fea)^^(a)^^(fe«) if) by 



aj Sij , 



(ka) 
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Note that with the coefficients given by Definition UJ applying the quantum Ito rule to 
Uf'^"'^ and uj:"\ defined in Definition [31 gives 




(9) 

to) 



Definition 5: Suppose that Assumptions [U [2l [3] and |4] hold. Let a be an element in C" 
and let i be an element in {1, . . . , n}. Define operators B^") and G\'^^ by 

^(") = A + -ajJ^Vy, 

=B + a,{W,, - + G^a, - oljG\W,,, 

g'^^^ =G, + a,W^J. 

Lemma 3: Suppose Assumptions [TJ [51 [3| and [HiioW. Let A, B, Y, Fi, Gi, Wij, K, Li and 
Sij for i,j e {1, . . . ,n} be the various operators occuring in Assumption [TJ [51 [3| and [^ 
Let K^°'^ and l\"^ fori €{!,..., n} be given by Definition\Mand let A("), and g|"^ 
for i e {1, . . . , n} be given by Definition\S[ Then 

-(G|"'-i^.yfiA("))Fo, (10a) 
if = Po (s'") - A(")yf ^A'")) Po, (10b) 

i.e. DeBnitionUl holds with A ^ A^°'\ B = B^"\ G, = G^"'^ L, = and K = if("). 
Moreover, Assumptions\^\^\^and\^hold for the altered coefficients with Pq and Y^^ 
unchanged. 

Proof. To show that Definition [T] holds for the altered coefficients, substitute G^"'' and 
^(") from Definition [5l and l\"-' from Definition [H into Eq. pOap . This gives 

L, + a.S^j = (l, + a,W^, + a,F,Yf^F^Wij,^ Pq, 

which holds if we substitute Sij = (^Wij + FiY^^ F^Wij^ Pq from Dcfinition[T] Further- 
more, substituting A^") and i?'"^ from Definition [5l and if*^"' from Definition [4] into Eq. 
(|10b|) gives 

ctiSijaj + CiiLi - ajL\Sij = PaaiWijajPa + P^GiPQai - ajP^GlWijPQ 
- Po{F,a, + A)yf i(A - a,F}Wi,)PQ + Po^i"f 'A/'o- 



13 



This holds if we can show that 



U = Po (G, - F.Yf^A) Po 

Lis,, = Po (gIw.j - AY-'F^W,,') Po. 



(11a) 
(lib) 
(11c) 



Equations (|llap and (jllbp are satisfied by Assumption 2] as PiLi — PiStj = 0. Note 
that Eq. (|llcp holds if we can show 



Now recall that Po{A + A'')Pi = -PoG^PjPi, and F + Ft ^ ~F^F, (see Eq. ©) by 
Assumptions m [3 and O Moreover, YY^-^ PiF^WijPo = PiF^Wi-jPo by Assumption [3] 
which shows that Eq. (jllcp is satisfied. 

We now show that Assumptions [1] [21 El and 21 hold for the altered coefficients, with Po 
and Y^^ unchanged. Assumption [H holds for the altered coefficients since, by Definition 
|3l we have t/j*^""* = W{ft)'^ul^^W{ft) which is clearly unitary. By Assumption |2l for the 
original coefficients and Definition|4land|5l we see that Assumption|2lholds for the altered 
coefficients. Assumption|3lon the altered coefficients is seen to hold by direct substitution 
of the coefficients in Definition IH and |5l followed by application of Assumption |3l for the 
original system. Assumption 01 holds if PiL^^^ = PiLi + aiPiSij — 0, which follows from 
Assumption |4l on the original system. □ 

Lemma |3l shows that Proposition [H holds with T^''"'^ and T^"'' replacing T^'^'' and Tt, 
respectively. 

Corollary 1: Suppose that Assumption [Tl|2l|3l and IH hold. Let a be an element of C". 
We have 



for aU X e B{no) and < s < oo. 

Proof of Theorem mH Let t > 0. Let / be a step function in L'^{[0,t\;C), i.e. there 
exists an m G N and ^ to < ti < . . . < tm — t and ai , . . . , e C" such that 



lI (^5u + F.Y-^f}^ WijPo = PoG\WijPo - PoAY-^f}Wi,Po. 



Substituting Li from Definition [1] this becomes 



PoA^Y-^^ f}WijPo + PoG\f,Y^^f}WijPo 



PoA^Y-^^f}f,Y^^f}WijPo+PoAY^^f}Wi,Po = 0. 




S £ [ti-i,ti) 



f{s) = ai, 



Vi e {1, . . . , m}. 
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The cocycle property of solutions to QSDE's and the exponential property of the sym- 
metric Fock space lead to 

r/'^) {X) = rif . . . T^^'ZL (^) ' ^ e B{n) , 

It is easy to see that Corollary [T] also holds for the difference of a finite product of maps 
T^^^^^_^_^ and a finite product of maps r(^°L']._^. This leads to 



lim 



lim 



Tt^''^\x) ^ Ty\x) 



rp(kam) 



T^'^:;^_^ix)-n^->...nx_^ix) 



-n(ai) 



= 0, 



X e B{Ho). 



This immediately yields for all step functions / G i^([0, t]; C") and v E Hq 

lim J/f (g) 7r(/) = Utv ® 7r(/). 



(12) 



Note that the step functions are dense in i^([0, t]; C"). This means that Eq. holds 
for aU / e L'^{[0, t];C"). Now note that for aU / e L'^{R+; C") and t < s < oo, we have 

(e.g. m) 

This means that the result in Eq. is true for aU / e L^(M+; C"). We now have 

lim [//''V = C^tV', 



for all ip in 9 = span{w «) 7r(/); v £ Hq, f e L2(R+;C")}. Theorem O then follows 
from the fact that is dense in Tip <Xi (e.g. (20]). 



5 Discussion 



In this article we have studied adiabatic elimination in the context of the quantum 
stochastic models introduced by Hudson and Parthasarathy. We have shown strong 
convergence of a quantum stochastic differential equation to its adiabatically eliminated 
counterpart, under four assumptions. Physically, the first Assumption [T] enforces the 
unitarity of the initial QSDE model. Assumption [2] ensures an appropriate scaling in 
the coupling parameter k such that we can distinguish excited and ground states in our 
system. Assumptions |3] and |4] ensure the existence of a limit dynamics independent of 
k. Note that Assumption [4] specifically forbids any quantum jumps which terminate in 
an excited state, the presence of which would preclude the construction of a valid limit 
dynamics. 

Although a Dyson series expansion for t/j'^-' (e.g. in terms of Maassen kernels [18]) would 
provide a lot of intuition for the results we have obtained (see [13] and [S] Chapter 5, 
Section 4]), we have chosen a proof along the lines of semigroups and their generators. 
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An infinitesimal treatment has the advantage that it can exploit the existence of results 
such as the quantum Ito rule J3] , the Trotter-Kato Theorem [^TJ [TB] and the technique 
due to Kurtz [T7j . 
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